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We report on a recent calculation of the two-loop virtual QCD correc- 
tions to the W boson pair production in the quark-anti-quark-annihilation 
channel in the limit where all kinematical invariants are large compared to 
the mass of the W boson. Our result is exact up to terms suppressed by 
powers of the W boson mass. The infrared pole structure is in agreement 
with the prediction of Catani's general formalism for the singularities of 
two-loop QCD amplitudes. 
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1. Introduction 

W pair production via quark-anti-quark-annihilation is a very important 
process at the Large Hadron Collider (LHC) for two main reasons. Firstly, 
it can serve as a signal process in the search for new Physics since it can 
be used to measure the vector boson trilinear couplings as predicted by the 
Standard Model (SM). Deviations from these predictions may come from 
decays of new heavy particles into SM vector boson pairs or anomalous 
couplings and would signal new Physics [1]. Secondly, qq — > W + W~ is the 
dominant irreducible background to the promising Higgs discovery channel 
pp — > H — > W*W* —> lul'v' in the mass range Mfn ggs between 140 and 180 
GeV [2j. In this paper, we describe the computation of the two-loop virtual 
QCD corrections to the W boson pair production in the quark-anti-quark- 
annihilation channel. The full result will appear in [3]. 

Due to its importance, the study of W pair production in hadronic 
collisions has attracted a lot of attention in the literature. The Born cross 
section was computed almost 30 year ago [4], whereas the next-to- leading 
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order (NLO) QCD corrections to the tree-level were computed in Refs. 
[?1 [U [9] and were proven to be large. They enhance the tree-level by 
almost 70% which falls to a (still) large 30% after imposing a jet veto. 
Another process that adds to the pp — » W + W~~ background is the W pair 
production in the loop induced gluon fusion channel, gg — > W + W~. This 
contributes at 0(a^) relative to the quark-anti-quark-annihilation channel 
but is nevertheless enhanced due to the large gluon flux at the LHC. The 
corrections from gluon fusion increase the W pair background estimate by 
almost 30% after certain experimental Higgs search cuts are imposed [10, 

im. 

Given that the NLO QCD corrections to the background are large and 
also that the cross section for the process H — > WW — > Ivl'v' (signal process 
for the Higgs discovery) is known at NNLO [12], the NNLO corrections to 
qq — > W + W~ need to be computed. This will allow a theoretical estimate 
for W production from gg-annihilation with accuracy better than 10%, as 
well as having both the signal and the background process calculated at the 
same order of the perturbative expansion (NNLO). In this paper, we address 
the task of computing the NNLO two-loop virtual part, more precisely the 
interference of the two-loop with the Born amplitude. We work in the limit 
of fixed scattering angle and high energy, where all kinematical invariants 
are large compared to the mass m of the W. 



2. The calculation 

Our methodology for obtaining the massive amplitude (massless fermion- 
boson scattering was studied in Ref. [13]) is very similar to the one followed 
in Refs. [HI [151 US], an evolution actually, of the methods employed in 
Refs. |17} 118] . The amplitude is reduced to a form that contains only a 
small number of integrals (master integrals) with the help of the Laporta 
algorithm [19] . In our calculation there are 71 master integrals. Next comes 
the construction, in a fully automatic way, of the Merlin-Barnes (MB) repre- 
sentations [2U\ I2T] of all the master integrals using the MBrepresentation 
package [22] . The representations are then analytically continued in the 
dimension of space-time with the help of the MB package [23J revealing the 
full singularity structure. An asymptotic expansion in the mass parameter 
is performed by closing contours and the integrals are finally resummed, 
either with the help of XSummer [23], or the PSLQ algorithm [25] . 

We will present here the result for a non-planar two-loop box master 
integral with two massive legs (Fig. [I}. The Mellin-Barnes representation 
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Fig. 1. A non-planar two-loop box master integral with two massive legs, pj are 
the external momenta, with p\ = p 2 = and p 3 = p 4 = m 2 , m the mass of the W. 



of this integral is 7-fold and reads: 
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This specific master integral is needed up to order 0((m? / s) -1 ). After the 
analytical continuation in the spacetime dimension and the expansion in the 
mass one is able to perform the resummation which yields: 



Inp — + 



1 



2m.se 4 



H 3 {~ 2L m + L x + L^} 



89vr 2 



4L? 



+ 



4L x L m 



4L y L m + L 2 + L 2 + 2L X L 



16IA 



12 ] m s e I 3 
~~ 4L 2 L m — 8La;LyL m + 



+ 8L x L m + 8L y L ri 



4L 2 L m 



897r 2 L m 



+ 



2li 2L 



+ -^ + 2L x L 2 



SOL^TT 2 SQLyTT 2 



+ 



1 16Lf 



32L a; L^ l 



4 



USTRON07 PRINTED ON FEBRUARY 2, 2008 



32L^L^ 2t 2 ot 2t 2 1 t t 2 1787T L m 
3 1" ° L x L m + ° L j/ L m + lbL x L y L m 3 

8L^.L m 8LyLm. 2 2 160C3Lm 



+ 



3 



178La;7r 2 L m 178L y 7r 2 L m L 4 L„ 4L X L^ 2 
+ § + § + y + y + -3- + 2L x L y 

+ 4L 3 x L y 80L X ( 3 80L y ( 3 IIIItt 4 89L 2 tt 2 



720 6 



89L 2 vr 2 SQL^tt 2 
6 3 



(2) 



We have used above the compact notation m s = m 2 /s, L m = log (m 2 /s), 
L x = log(x) and L y = log(l — x). s = (pi + P2) 2 and t = (pi — P3) 2 are 
the usual Mandelstam variables and s = t/s. One of the checks in our 
calculation was to verify numerically the analytic results for the masters. 

The renormalization of the amplitude involves only the renormalization 
of the strong coupling constant a s \26\ I27j . A non-trivial check for our 
computation is that the infrared pole structure of our renormalized result 
is in agreement with the prediction of Catani's general formalism for the 
singularities of the two-loop QCD amplitudes as described in Ref. [28] . 



3. Outlook 

We have calculated the two-loop virtual QCD corrections to the W boson 
pair production in the quark-anti-quark-annihilation channel in the high 
energy limit. The two- loop result, along with the square of the one-loop 
2^2 process have to be combined with the tree-level 2^4 and the one- 
loop 2 — > 3 processes in order to obtain physical cross sections. Combining 
all these contributions will enable the analytic cancellation of the remaining 
infrared divergences. Initial state singularities will have to be absorbed into 
parton distribution functions of the hadrons (protons) in order to match 
with a precise parton evolution at NNLO [29j [30] . 
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